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We describe aspects of particle creation in strong fields using a 
quantum kinetic equation with a relaxation-time approxima- 
tion to the collision term. The strong electric background field 
is determined by solving Maxwell's equation in tandem with 
the Vlasov equation. Plasma oscillations appear as a result of 
feedback between the background field and the field generated' 
by the particles produced. The plasma frequency depends 
on the strength of the initial background field and the colli- 
sion frequency, and is sensitive to the necessary momentum- 
dependence of dressed-parton masses 
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Ultra-relativistic heavy-ion collisions are complicated processes and their 
understanding requires a microscopic modelling of all stages: the formation, 
evolution and hadronisation of a strongly coupled plasma. If the energy density 
produced in the interaction volume is large enough, then the relevant degrees of 
freedom are quarks and gluons. The terrestrial recreation of this quark gluon 
plasma (QGP) will aid in understanding phenomena such as the big bang and 
compact stars. 

Construction of the Relativistic Heavy Ion Collider at the Brookhaven Na- 
tional Laboratory is complete and the initial energy density: e ~ 10 — 100 
GeV/fm 3 , expected to be produced in the collisions at this facility is certainly 
sufficient for QGP formation. Experimentally there are two parameters that 
control the conditions produced: the beam/target properties and the impact 
parameter. Varying these parameters changes the nature of the debris mea- 
sured in the detectors. Signals of QGP formation and information about its 
detailed properties are buried in that debris jl). Predicting the signals and 
properties requires a microscopic understanding of the collisions, including their 
non-equilibrium aspects. 

In the space-time evolution of a relativistic heavy ion collision the initial 

*To appear in the Proceedings of Quark Matter in Astro- and Particlephysics, a workshop 
at the University of Rostock, Germany, November 27 - 29, 2000. Eds. D. Blaschke, G. Burau, 
S.M. Schmidt. 



state is a system far from equilibrium. This system then evolves to form an 
equilibrated QGP, and the investigation of that evolution and the signals that 
characterise the process are an important contemporary aspect of QGP research. 

The formation of a QGP is commonly described by two distinct mechanisms: 
the perturbative parton picture 0] and the string picture || . In the parton pic- 
ture the colliding nuclei are visualised as clouds of partons and the plasma 
properties are generated by rapid, multiple, short-range parton-parton interac- 
tions. In the string picture the nuclei are imagined to pass through one another 
and stretch a flux tube between them as they separate, which decays via a 
nonperturbative particle-antiparticle production process. These approaches are 
complimentary and both have merits and limitations. Once the particles are 
produced the subsequent analysis proceeds using Monte-Carlo event generators 

[H§. 

Herein we employ the nonperturbative flux tube picture M. A flux tube is 
characterised by a linearly rising, confining quark- antiquark potential: V q q{r) = 
a r. The string tension can be estimated in lattice simulations using static 
quark sources, which yields a ~ 4Aq CD ~ 1 GeV/fm. This string tension can 
be viewed as a strong background field that destabilises the vacuum and the 
instability is corrected through particle-antiparticle production via a process 
akin to the Schwinger mechanism ||. Figure |l| is an artist's impression of this 
process. 

Assuming a constant, uniform, Abelian field, E, one is able to derive an ex- 
pression for the rate of particle production via this nonperturbative mechanism 



cxp 



27r(m 2 +p]_) 



(1) 



where m is the mass and e the charge of the particles produced. It is plain from 
this equation that the production rate is enhanced with increasing electric field 
and suppressed for a large mass and/or transverse momentum. 

The production of charged particles leads naturally to an internal current. 
That current produces an electric field that increasingly screens and finally com- 
pletely neutralises the background field so that particle creation stops. However, 
the current persists and the field associated with that internal current restarts 
the production process but now the field produced acts to retard and finally 
eliminate the current. . . . This is the back-reaction phenomenon and the natu- 
ral consequences are time dependent fields and currents || . One observable and 
necessary consequence is plasma oscillations. Their properties, such as frequency 
and amplitude, depend on the initial strength of the background field and the 
frequency of interactions between the partons. It is clear that very frequent 
collisions will rapidly damp plasma oscillations and equilibrate the system. 

The process we have described can be characterised by four distinct time- 
scales: 

I. the quantum time, r qu , which is set by the Compton wavelength of the 
particles produced and identifies the time-domain over which they can be 
localised/identified as "particles;" 
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Figure 1: Left: For E = the vacuum is characterised by a completely filled 
negative-energy Dirac sea and an unoccupied positive-energy continuum, sepa- 
rated by a gap: 2st = 2 (m 2 +p 2 L ) 1 / 2 . Right: Introducing a constant external 
field: E = e x E, which is produced by a potential: A = —Ex, tilts the energy 
levels. In this case a particle in the negative-energy sea will tunnel through 
the gap with a probability ~ exp(— ire^/eE). Succeeding, it will be accelerated 
by the field in the — x-direction, while the hole it leaves behind will be accel- 
erated in the opposite direction. The energy-level distortion is increased with 
increasing E and hence so is the tunneling probability. eE can be related to the 
flux-tube string-tension. 



2. the tunneling time, r tM , which is inversely proportional to the flux tube 
field strength and describes the time between successive tunnelling events; 

3. the plasma oscillation period, t p i, which is also inversely proportional to 
field strength in the flux tube but is affected by other mechanisms as well; 

4. and the collision period, t co u, which is the mean time between two partonic 
collision events. 

Each of the time-scales is important and the behaviour of the plasma depends 
on their relation to each other pc|jll[ ; e.g., non-Markovian (time-nonlocality) 
effects are dramatic if r qu ~ r tu . Herein we focus on the interplay of the larger 
time scales and consider a system for which r p ; ~ t co u . 

Particle creation is a natural outcome when solving QED in the presence of 
a strong external field and a formulation of this problem in terms of a kinetic 
equation is useful since, e.g., transport properties are easy to explore. The 
precise connection between this quantum kinetic equation and the mean field 
approximation in non-equilibrium quantum field theory p2| is not trivial and 
the derivation yields a kinetic equation that, importantly, is non-Markovian in 
character ||,||J||@. 



The single particle distribution function is defined as the vacuum expectation 
value, in a time-dependent basis, of creation and annihilation operators for single 
particle states at time t with three- momentum p: at(t), ap(t); i.e., 

ffit):={0\aUt)ap(t)\0). (2) 
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The evolution of this distribution function is described by the following quantum 
Vlasov equation: 



df±(P,t) df ± (P,t) „ M df ± (P,t) 

i i* 

= =W±(t) / dt'W±(t') x [l±2f±(P,t')]co S [x(t',t)} + C ± (P,t), 

* •* — oo 

where the lower [upper] sign corresponds to fermion [boson] pair creation, C 
is a collision term and W± are the transition amplitudes. The momentum is 
defined as P = (pi,P2, P\\ (t)), with the longitudinal [kinetic] momentum P|| (t) = 
p 3 -eA(t). 

We approximate the collision-induced background field by an external, time- 
dependent, spatially homogeneous vector potential: A^, in Coulomb gauge: 
A n = 0, taken to define the z-axis: A — (0, 0, A(t)). The corresponding electric 
field, also along the z-axis, is 

E(t) = -A(t) = . (4) 
For fermions [ p^Hf5[ ] and bosons Jn],|l5) the transition amplitudes are 



where the transverse energy e± = y/ m 2 + p\ , p± = (p%,P2), and oj(t) 
e\ + PfAt) is the total energy. In Eq. (|), 



x(t',t) = 2[Q{t)-Q{t')], 9(t)= / dt'uj(t'), (6) 



is the dynamical phase difference. 

The time dependence of the electric field is obtained by solving the Maxwell 
equation 

- A ± (t) = £±(t) = -j°*(t) - j cond (t) - j pol (t) (7) 

where the three components of the current are: the external current generated, 
obviously, by the external field; the conduction current 

3Wt ) = P ± e/^J|/ ± (P,*), (8, 

associated with the collective motion of the charged particles; and the polarisa- 
tion current 



j P oi(t) = g±t 



d 3 P P\\(t) 
(2tt) 3 w (p )f ) 



S(P,t) eE^P^t) 



W±(P,t) 8^(P,t) 



, (9) 
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Figure 2: Time evolution of the fermion's electric field obtained with three 
different initial field strengths. The amplitude and frequency of the plasma 
oscillations increases with an increase in the strength of the external impulse 
field, Eq. (|l0|). The reference time-scale is the lifetime b of the impulse current. 
(The field and time are given in arbitrary units, and b = 0.5.) 



which is proportional to the production rate. Here <?_ = 2, g + — 1 and all fields 
and charges are understood to be fully renormalised, which has a particular 
impact on the polarisation current fl7| . 

We mimic a relativistic heavy ion collision by using an impulse profile for 
the external field 

E ex = -^-sech 2 (t/b), (fO) 

which is our two-parameter model input: the width b and the amplitude A 
are chosen so that the initial conditions are comparable to typical/anticipated 
experimental values. This is the seed in a solution of the coupled system of 
Eqs. (|) and (0)-@, and in Fig. || we illustrate the result for the electric field 
in the absence of collisions. 

The qualitative features are obvious and easy to understand. The external 
impulse (the collision) is evident: the portion of the curve roughly symmetric 
about t = 0, where the field assumes its maximum value A$/b. It produces 
charged particles and accelerates them, producing a positive current and an 
associated field that continues to oppose the external field until the net field 
vanishes. At that time particle production ceases and the current reaches a 
maximum value. However, the external field is dying away: it's lifetime is t ~ 
0.5, so that the part of the electric field due to the particles' own motion quickly 
finds itself too strong. The excess of field strength begins to produce particles. 
It accelerates these in the opposite direction to the particles generating the 
existing current whilst simultaneously decelerating the particles in that current. 
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That continues until the particle current vanishes, at which point the net field 
has acquired its largest negative value. Particle production continues and with 
that a negative net current appears and grows. . . . Now a pattern akin to 
that of an undamped harmonic oscillator has appeared. In the absence of other 
effects, such as collisions, it continues in a steady state with the magnitude and 
period of the plasma oscillations determined by the two model parameters that 
characterise the collision. 

The quantum kinetic equation, Eq. (||) , describes a system far from equilib- 
rium and therefore any realistic collision term C valid shortly after the impact 
must be expected to have a very complex form. However, as Fig. |2| illustrates, 
the evolution at not-so-much later times is determined by the properties of the 
particles produced and not by the violent nonequilibrium effects of the collision. 
This observation suggests that the parton plasma can be treated as a quasi- 
equilibrium system and that the effects of collisions can be represented via a 
relaxation time approximation [[l7|Jl^,p^,po|j2l|| : 



C(p,t) 



-(*) 
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(ii) 



with t{£) an in general time-dependent "relaxation time" (although we use a 
constant value r(t) = r r in the calculations reported herein) where 



fJ(p,T(t)) = 



exp 



p v W{t) 
T(t) 



1 



(12) 



is the quasi-equilibrium distribution function for fermions. Here T(t) is a local- 



temperature and u v (t), 



1, is a hydrodynamical velocity 



(With our 



geometry, u v {t) = (1, 0, 0, v{t))/[l - v 2 ^)] 1 / 2 .) 

Our definition of quasi-equilibrium is to require that at each t the energy 
and momentum density in the evolving plasma are the same as those in an 
equilibrated plasma; i.e., we require that 



e f (t) = e^(t),p f (t)=^(t) 
where, as one would expect, 



'(*) = 



d 3 p 



d 3 p 



(13) 



p(t)f- q (P,t), (14) 



and 



where 



d 3 p 

d 3 p 
(2^)3 



U(j>, t) [/_ (p, t) - Zf_ (p, t)] 

p(t)[f-(p,t)~zf-(p,t)] , 



( e£± 
\4uj 3 



E{tf 



o-2t/r r 



dt'E{t') 2 e 2t 'l^ 



(15) 
(16) 

(17) 
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Figure 3: Time evolution of the electric field for three different relaxation times. 
As one would anticipate, collisions damp the plasma oscillations: their frequency 
and amplitude decreases with decreasing relaxation time. Solution obtained 
using the impulse profile, Eq. (|f0|), with Aq = fO, b = 0.5. The field and the 
time are given in arbitrary units. 



is a regularising counterterm. Adding Eqs. (|13|) to the system of coupled equa- 
tions embeds implicit equations for the temperature profile, T(t), and collective 
velocity, 

Solving the complete set of coupled equations is a straightforward but time 
consuming exercise. For the present illustration we employ the minor simplifi- 
cation of assuming that the equilibrium energy density is that of a two-flavour, 
massless, free-quark gas; i.e., 

t eq (T(t)) = ^T 4 (t) (18) 

and then proceed. 

The solution obtained for the electric field using the impulse profile, Eq. (p"o|), 
is depicted Fig. ||. Here the behaviour is analogous to that of a damped oscil- 
lator: the frequency and amplitude of the plasma oscillations diminishes with 
increasing collision frequency, l/r r (friction). 

In the examples presented hitherto we have employed a constant fermion 
mass. However, in QCD the dressed-quark mass is momentum-dependent p2| , p3[ 
and that momentum-dependence is significant when mo < Aqce>, where m is 
the current-quark mass. That can be illustrated using a simple, instantaneous 
Dyson-Schwinger equation model of QCD, introduced in Ref. |Q, which yields 
the following pair of coupled equations for the scalar functions in the dressed- 

A shortcoming of the approach presented thus far is that it neglects the possibility of dis- 
sipative inelastic scattering events transforming electric field energy directly into temperature. 
However, we have almost completed an extension valid in that case. 
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Figure 4: The scalar functions characterising a dressed-u-quark propagator, as 
function of momentum, obtained using the simple DSE model introduced in 
Ref. |§, see Eqs. |L9|), @. We plot the functions as obtained at T = T c = 
0.17 GeV, which is the critical temperature for deconfinement in the model fl25{. 



quark propagator: S(p) = l/[«7 ■ pA(p 2 ) + B(p 2 )], 

mt) - m "^7wwBwm li - 21 - m ' h (19) 

MP,t) = -™-y (20) 
m + B(p,t) 

where rj is the model's mass scale. In this preliminary, illustrative calculation 
we discard the distribution function in Eq. (|l9|). 

The solution obtained using a current-quark mass mo = 5 MeV and with 
r) = 1.33 GeV is depicted in Fig. This value of the mass-scale parameter 
can be compared with the potential energy in a QCD string at the confinement 
distance, V qq {r = 1 fm) = ar ~ (2A QC d) 2 (1/Aqcd) = 4A QCD - 1.0 GeV. The 
dressed-quark mass function is m(p 2 ,T) = B(p 2 ,T)/A(p 2 ,T), and m(0,T) 
provides a practical estimate of the T-dependent constituent-quark mass |p6f : 
in this example m(0, T c ) sa 0.35 GeV. 

In our flux tube model for particle production we produce fermions with 
different momenta and following this discussion it is clear that the effective 
mass of the particles produced must be different for each momentum. That will 
affect the production and evolution of the plasma. To illustrate that we have 
repeated the last calculation using m(p) wherever the particle mass appears in 
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Figure 5: Electric field as function of time for Ao = 10.0, t = 5.0, b = 0.5. We 
compare the results obtained using particles produced with a constant mass, 
m = 1, with those obtained when the mass is momentum-dependent. For il- 
lustrative simplicity, we use the profile in Fig. |J but normalised such that, 
m(0) = m. The mass-scale is arbitrary. 

the system of coupled equations for the single particle distribution function. 
The effect on the electric field is depicted in Fig. 0, wherein it is evident that 
the plasma oscillation frequency is increased when the momentum-dependent 
mass is used because many of the particles produced are now lighter than the 
reference mass and hence respond more quickly to changes in the electric field. 

We have sketched a quantum Vlasov equation approach to the study of par- 
ticle creation in strong fields. Using a simple model to mimic a relativistic heavy 
ion collision, we solved for the single particle distribution function and the as- 
sociated particle currents and electric fields. Plasma oscillations are a necessary 
feature of all such studies. We illustrated that the oscillation frequency depends 
on: the initial energy density reached in the collision, it increases with increas- 
ing energy density; and the particle-particle collision probability, it decreases as 
this probability increases; and that it is sensitive to the necessary momentum- 
dependence of dressed-particle mass functions. The response in each case is 
easy to understand intuitively and that is a strength of the Vlasov equation 
approach. 
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